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Abstract 



A full selfconsistent set of equations is deduced to describe the kinetics 
and dynamics of charged quasiparticles (electrons, holes etc.) with arbitrary 
dispersion law in crystalline solids subjected to time-varying deformations. 
The set proposed unifies the nonlinear elasticity theory equation, a kinetic 
equation for quasiparticle excitations and Maxwell's equations supplemented 
by the constitute relations. The kinetic equation used |]l|j2| is valid for the 
whole Brillouin zone. It is compatible with the requirement for periodicity 
in /c-space and contains an essential new term compared to the traditional 
form of the Boltzmann equation. The theory is exact in the frame of the 
quasiparticle approach and can be applied to metals, semiconductors, as well 
as to other crystalline solids including quantum crystals and low-dimensional 
lattice structures. 
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I. INTRODUCTION 



There are two fundamental problems when dealing with quasiparticles in crystalline 
structures. The first one is related to the fact that the quasimomentum k and the dispersion 
law e(k) of a quasiparticle are well defined only in an ideal periodic lattice. In such a lattice 
the dispersion law as well as all other physical quantities are periodic functions in the 
reciprocal space (/c-space). However, in any real system the crystal lattice is deformed (e.g. 
by impurities, elastic deformations, external fields etc.). The most complicated problems 
concern the quasiparticle dynamics in a time-varying deformed crystal lattice. In an exact 
description all physical characteristics of a quasiparticle have to be periodic functions of the 
quasimomentum with periods which are functions of the coordinates and the time. This 
leads to a dependence of the Brillouin zone boundaries not only on the deformation at a 
given instant, but also on the local lattice velocity 

As for the stationary (or quasistationary) case, this difficulty has usually been passed 
over by introducing a local dispersion law e(k,Uik) and further expansion in powers of the 
small deformation tensor components Uik'- 



and u is the deformation vector. Such an approach (known as the local lattice approach) 
is based on the assumption that the deformations are small and smooth functions of po- 
sition and time and possesses all shortcomings of any linearized theory. In addition, the 
kinetic equation becomes incompatible with the requirement of the periodicity. To avoid 
any misunderstandings, note that Eqs.(l) and (2) are written in the co-moving frame. 

The second problem is related to the fact that the crystal lattice plays the role of a 
privileged coordinate frame and no Galilean transformations for quasiparticle characteris- 
tics exist. The lack of transformation laws for such quantities as energy, Hamiltonian and 




(1) 



where Xik are the deformation potential components. 




(2) 
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quasimomentum means, in fact, that there is no consistent quasiparticle mechanics. The 
most fundamental quantity, the dispersion law, is known in a co-moving frame attached to 
the lattice and this frame is even not inertial at time-varying deformations. On the other 
hand, all fundamental physical equations, such as conservation laws, variational principles, 
kinetic eqs. etc. take place in the laboratory frame (L-system). However the concept of dis- 
persion law does not exist in L-system. Hence, in principle, even if the mechanics equations 
for a quasiparticle were known in C-system, they remain unknown in L-system. 

These problems are so old as the quasiparticle approach itself is. They are well known, for 
example, in the theory of metals where many attempts have been made to derive a complete 
system of dynamic equations, consisting of equations from the theory of elasticity, a kinetic 
equation for the electron gas and Maxwell's equations (cf. Refs. and the bibliography 
cited there). 

The attempts to deduce the equations of motion for a quasiparticle in L-system can 
be divided into two groups corresponding to the twofold role of the dispersion law. In the 
co-moving frame it coincides with both energy and Hamiltonian. Therefore transformations 
typical of energy and Hamilton function have been proved. 

When considering e(k, Uik) as Hamiltonian a transformation by a substitutional function 
of the form 



where r = r' + u are the coordinates in L-system and p is supposed to be the corresponding 
quasimomentum. However, $ does not depend on the bare mass, m, of the quasiparticle 
and therefore does not take into account any inertial effects (e.g. the Stewart-Tolman effect 
in metals, centrifugal forces in rotating bodies etc). 

If £:(k, Uik) is considered as quasiparticle energy then the Galilean transformation applies 




(3) 



has been used and the following relations have been obtained (see e.g. 0,^): 



k = p + ^(up), H{p, r, t) = e(k, r', t) + up. 
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as a consequence of the requirement that energy density / ed^k (as a macroscopic quantity) 
has to obey Gahlean transformations. This yields 

8 = eo(k) + \ik{k)uik + (4) 

It has been shown by many authors that such a transformation is incompatible with the 

Boltzmann equation. That is why some artificial ways have been used the most successful 

being that of Landau (cf the footnote in Ref. 0). He has suggested that (in order to take 

into account the noninertial properties of the lattice frame) one has first to add to the 

. de 

dispersion law (1) the term — mu— , setting 

ak 

de 

e{k, Uik) = eo{k) + \ik(k)uik - '^^g^- (5) 

and then to apply the transformation by the substitutional function (3). The result is: 

H{p, r, t) = eo(k) + Xikik^Uik + - ^ (6) 

This procedure has been used in the most of recent works. However it cannot be well- 
grounded due to its internal inconsistency. In fact, if we consider the expression (5) as energy, 
and take a constant velocity u then we come to the wrong conclusion that the energy in an 
inertial frame could depend on the frame velocity. The same confusion follows for the energy 
density which must strictly obey the Galilean principle. Note, that we may not consider (3) 

as a Hamiltonian, because the Hamilton function in C-system coincides with the dispersion 

de 

law in accordance with the Hamilton equation r' = v = — — . 

OK 

These shortcomings have been removed in our previous works |l],0J3l (^^^ [00) where 
transformations to replace the Galilean ones have been deduced and the quasiparticle me- 
chanics equations in Hamiltonian form have been presented. This enabled us to derive a 
Boltzmann-like kinetic equation valid in the whole Brillouin zone. We have deduced a self- 
consistent set of equations for electrons in metals, taking into account some special features 
of the problem as the quasineutrality condition and the neglection of the displacement cur- 
rent in Maxwell's eqs. The magnetic permeability /i has also been taken constant (/i = 1). 



These approximations were good enough to develop the electron plasma hydrodynamics in 
crystalline metals as well as to consider magnetohydrodynamic effects [Q. 

The problem is more complicated in bad conductors and semiconductors, as well as at 
higher frequences, when the displacement current cannot be neglected, the quasineutrality 
condition does not hold and e and /i are functions of the deformation. We give in this work 
the solution of this general case. The only assumption is that the constitute relations are 
linear, i.e. A = hkEk, Bi = jliuHk- 

In the present work we shall consider electrons having in mind that the theory is valid for 
arbitrary quasiparticles. We suppose for simplicity that the crystal considered is isotropic 
in its undeformed state. This means that e and are taken scalars, but their derivatives 
with respect to coordinates are matrices and depend on the lattice symmetry. It is easy 
to generalize all the results for the case where e and /i are tensors of second rank. No 
essentially new results should be obtained for this case, but the corresponding relations are 
more cumbersome. 

Finally, we would like to note that the problem considered is related to some old questions 
about the electromagnetic forces acting to the body, some specific features of electrodynamics 
in moving media, the form of electromagnetic stress tensor in condensed matter, the role of 
momentum and quasimomentum etc. Since we derive the full set of dynamics equations all 
the forces are taken into account in a selfconsistent way. 

The paper is organized in the following way. In Sec. 2. we introduce new variables 
(discrete coordinates and ivariant quasimomentum), write evolution equations for lattice 
vectors in the real and reciprocal spaces and express in new notation metrical tensors, 
deformation tensor etc. In Sec. 3. we reproduce briefly some of our previous results on 
the Hamilton eqs. and kinetic equation we need for this work. In Sec. 4 we deduce the 
full set of equations which describes the behaviour of quasiparticles in deformable solids in 
electromagnetic fields. 
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II. NOTATION 



Following our previous works |]T],0 (see also we consider a lattice with primitive 

vectors Sia, a = 1, 2, 3 and introduce discrete coordinates N"" as the number of steps (each 
being equal to slo) in the lattice from the origin to a given point. In this notation the differ- 
ential coordinates dr which are considered as physically infinitesimal (i.e. large compared 
to the lattice constants but small compared to any macroscopic distance of interest) can be 
written in the form: 

dr = a„dA^° + udt 

or 

dN"" = a"c/r - a^iidt 
where a" are the reciprocal lattice vectors which satisfy the relations 

aQ.a'^ = 5^, CLaiO-k ~ ^ik- i'^) 

It is seen from the above equations that 

3r 

a- = ViV°, ^o.= g^, u=-a„iV". (8) 

The time-evolution equations for the vectors a^ and a" can be deduced from plain geomet- 
rical considerations (Appendix 1) and written as follows 0: 

a^, + (uV)a^ - (a„V)u = (9) 

a" + V(a°u) = 0. (10) 
In the notation used the metrical tensors in the real and reciprocal space are, respectively: 

ga/3 = a^a^, = sl^sl^ (11) 

Then the lattice cell volume equals g^^'^, where g = det ga/s- 

The relations between the components of the metrical tensors and the deformation tensor 
Uik follow from the expression for the interval: 

6 



ds'^ = QapdN^'dN^. (12) 
The squared interval between the same points in the undeformed crystal is 

Hence 

ds^ - dsl = 2wo,pdmdN'^ (13) 

where 

1 

Wap = ^{9al3- dap) (14) 

plays the role of the deformation tensor in our notation. The invariance of the interval 
yields: 

WapdN'^dN^ = Wikdx'dxK (15) 

Taking into account relations (7) and (8) one obtains 

Wik = Wapa'-al- (16) 

In order to obtain the components Uik of the tensor of small deformations (2) as well as 
to find the relation between the deformation vector u = r — ro and the discrete coordinates 
N"" introduced, let us note that the quantities dui/dx^ obviously coincide with the matrix 
elements which describe the coordinate transformations 

Xi^Xi^aikXk (17) 

and, consequently, the lattice vector transformations: 

o o t io\ 

dai —O'ai -\-OCik Cbak ■ (18) 

If the deformations are small, then, 

aa = Sia +S3La, = + 5a^. (19) 



Multiplying these two equations and taking into account relations (7) and (17) yield in 
a linear approximation with respect to 5a.^ 

^ -{aj3i 5al) a^k^ aik a^k ■ (20) 

On the other hand the change of the discrete coordinates owing to the deformation can 
be written in the form: 

AT" = 7V« - w". (21) 

where is the deviation from the value Nq in the ideal undeformed lattice. Taking into 
account (8) one obtains 

a« ^ VA^" = a" - Vw"" (22) 

and hence, 

= -IS-- (^^' 

It follows from (20) and (23) that 

dui o dw"^ d Clai , , 

- — = aik = o.ai — = — ^ (24) 

OXk OXk OXk 

and therefore the deformation vector u and the deformation tensor un^ in this linear case 

Till 1 1 1 ° 1 / "^Wq, dwR \ , 

are related to their discrete coordinate analogies, w and w^^^ - I + I > m the 

following way: 

u = a^, w , Uik= Wa(3 < (25) 

o 

where Wa = 9a/3 "U^ are the covariant components of w. 

In the same way one can obtain the full deformation tensor components [5] : 

1 / o 9u o 9u du du \ 
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III. DYNAMICS AND KINETICS OF QUASIPARTICLES 



The starting point when deriving Hamilton equations for quasiparticles is that in the 
co-moving frame (C-system) the dispersion law £{\i,g°'^) coincides both with Hamiltonian 
and with energy. Hence, the equations of motion in C-system can be written in the form 

., ds ■ de 

We have to determine Hamiltonian i?(p, r, t) as a function of the coordinates r and 
quasimomentum p in L-system in a way to have canonical equations: 

. dH . dH 



r 



(28) 



9p ' dr ' 

According to the general theory of Hamilton mechanics, the Hamiltonian and the mo- 
mentum can be obtained as derivatives of the action 5'(r, t) with respect to the time t and 
the coordinates r. Equations (27) show that the variables k and r' are canonically conju- 
gate. However, if one considers r' as a continuos variable, then the corresponding quantity k 
has to be considered as momentum. Quasimomentum has to be conjugate to some discrete 
coordinate. Such coordinates, N (with components N"')., were introduced in Sec. 2. Hence, 
we are able to take the action 5'(N, t) as a function of these coordinates and the time. Then 
we can determine a Hamiltonian H{k, N, t) and a new quasimomentum k as follows: 

We call K the invariant quasimomentum, because all physical quantities, written as functions 
of K, are periodic with a constant period 27r (not 27ra°). 

Let us consider the local dispersion law e(K, (yf"^) as a function of the invariant quasimo- 
mentum and the metrical tensor g"^ . Since it coincides with the Hamiltonian, one has 

Now let us consider a real electron which is executing a quasiclassical motion. Its wave 
function in the new variables has the form 
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V'(N,t)~exp|^5o(N,t)| (31) 

where 5*0 (N, t) is the classical action. The transformation law for the action follows from the 
transformation properties of the phase of the wave function under Galilean transformations 
[9]: 

S = So + mur — (32) 

where m is the mass of a free particle. 

The Hamiltonian and the quasimomentum in the Laboratory frame can now be obtained 
as follows 0: 

p=(^) =SL^ka + mu (33) 



t 



H{p,r,t) = -i^—j =e + pu-^ (34) 

where e = e{sLa{p — mii), (7°^) is a periodic function of p with periods 27r/ia" determined by 
the reciprocal lattice vectors corresponding to the deformed local lattice. This is the reason 
to call p the quasimomentum of the quasiparticle in L- system. 
It follows from (33) that 

ka = aa(p — mil) = ka^. (35) 

Hence, the invariant quasimomentum components are equal to the scalar product of the 
usual quasimomentum k in C-system and the primitive vectors of the locally deformed 
lattice. 

The energy £^ of a quasiparticle in L-system obeys the Galilean law: 

mu'^ . de mu^ . , . 

S = —— + mn— + € = —— + pou + e (36) 
2 ap 2 

de 

where po = rn—— is the average momentum (the mass flow) in C-system. 
ap 
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Equations (33)-(36) replace the Galilean transformations, which are not valid for quasi- 
particles because of the priviged of the crystal lattice frame. 
We are able now to write the Boltzmann equation: 

9/ ^ 9/ (37) 

dt dr dp dp \ dr J 

where / is the collision operator. Note, that this equation becomes well defined for quasi- 
particles only after obtaining the Hamiltonian and the Hamilton equations. 

It has been shown in |0 that this equation is compatible with the requirement of period- 
icity. This can be seen also from its form if /(k, r, t) is taken as a function of the invariant 
quasimomentum k and the quantities p and H are substituted from (33) and (34). This 
yields: 

d/ de fdf\ df \ du lde\ de , • fx /oon 

Tt^^-d^X-di)^-^-^\^Tt^\d.)^' "^^^^ ^ ^^^^ u - f| = // (38) 

d d 

where ^ ~ ^ + ('^^) ^^"^ ^ quantities are differentiated with respect to the coordinates 

and the time at constant k. 

The term m — - takes into account noninertial properties of the local frame. This is the 
dt 

term which is responsible for the Stewart-Tolman effect in metals. 
The term 

de 

m— — afl X curlii (39) 

dKfi 

is of essentially new kind and cannot be obtained in linear theories. It is proportional to the 
bare mass m and, hence, is also responsible for noninertial effects. 
The Lorentz force in our notation has the form 

p de p 
F = -eE - - x B - -u x B (40) 

C dHa C 

where E and B are the strengths of the electric field and magnetic induction, respectively. 
Substituting (40) into (38) yields the kinetic equation for charged quasiparticles with a 
charge — e ( e > 0): 
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d/^^ de _^ df \fde\ ^^^^^ ^ x b1 - // (41) 

where 

~ c dii ~ Tnc 

eE = eE+ -u X B + m— , B = B curlu. (42) 

c dt e 

The first expression in (42) consists of two parts — the electrical force eE' where E' — 

E H — u X B is the field in the co-moving frame, and the incrtial force m— . The second 
c dt 

~ c 

relation in (42) can also be written in the form B = curl (A mix) where A is the vector 

potential in an agreement with the well-known rule of replacement of the particle momentum 

g 

in a magnetic field V ^ V A where now V = mix . 

c 

We shall use also another form of the kinetic equation. The reason is that further on we 
need to integrate some physical quantities over the Brillouin zone, take integrals by parts 
as well as differentiate with respect to the coordinates and the time. However, the Brillouin 
zone boundaries are functions not only of the deformation at given instant, but also on the 
local lattice velocity. The integration over the Brillouin zone does not commutate with the 
differentiation with respect to the time and the coordinates. As a result some fluxes appear 
through the zone boundaries. This effect is important for nonequilibrium systems, open 
Fermi-surfaces as well as for other cases when the partition function or its derivatives do 
not vanish on the zone boundaries. This inconvenience can be passed over by introducing a 
renormalized partition function 

^{K,v,t)^f/^. (43) 
The kinetic equation for (^(k., r, t) has the form 

+ div I (u + a.^) ^} - a. A 1^ (^vs + m^) 

de 1 
— mio— — a/3 X curlii — Yip > = Ild. (44) 

J 

In this notation the differentiation with respect to t and r is carried out at constant 
K and hence commutates with / dP'n. Results obtained by such a procedure can easily be 
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rewritten in the previously adopted variables by the following substitutions 

d d d 



' dKa c?k dp 



<^(K,r,t)... = ((<^...)), (45) 



{2'Kh) 

d^r...= I d^N'^g^/'^ . . . (46) 



IV. CONSERVATION LAWS AND DYNAMICS EQUATIONS 

In order to avoid cumbersome expressions we shall make our consideration in three steps. 
First we write the conservation laws and consider the problem in a general form without 
taking into account the explicit form of the Maxwell's equations and electromagnetic forces. 
Then we consider Maxwell's equations and constitute relations for moving media and finally 
we combine the results and obtain the full selfconsistent set of equations. 

We start with the conservation laws. 

The continuity equation for quasiparticles is 

mh + div jo = (47) 

where 

—f\=m(—fj+mnu. (48) 

This equation follows directly from the kinetic equation 
The total mass current is 

Jo = pu + jo, (49) 

where 
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= (50) 

is the quasiparticle mass current with respect to the lattice and p = po + ^'t- is the full mass 
density written as a sum of the lattice mass density po and the quasiparticle mass density. 
The quantities p and Jq satisfy the mass continuity equation 

p + divJo = 0. (51) 

The full momentum J is a sum of Jq and the field momentum g: 

J = Jo + g. (52) 

Note, that in this case the full momentum does not coincide with the mass flow! 

Our aim is to determine momentum and energy fluxes Iljfe and Q in such a way as to 
satisfy the continuity equation (51), the momentum conservation law 

Ji + VfeHife = 0, (53) 

and the energy conservation law 

£; + div Q = 0. (54) 

The energy in L-system is given by the expression 

E = -^pou' + E,{g^P) + {{£^)) + W, (55) 

where Eo^g^^^) is the strain energy in C-system, and W is the field energy. 
The time derivative of the energy (55) is then 

E^pun+ \pu- + I {{s^)) + mua. ((^^)) 

The time derivative of the elastic energy Eq^q"^) can be taken in the following way. 

Eo = §^9"' = -aapa^a^ (57) 
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where 



^o^^ - -^d^- (58) 



Substituting af from the evolution equation (10) in (57) yields 



On the other hand 



^. = ...af(g«. + <g. (59) 



V.^. - ^V,- ^ (60, 



It follows from (59) and (60) that 

Eo = (ra(3at4^^-uVEo. (61) 
The time derivative of the full momentum (52) gives 

= - j + pu + pu + ma. ((^^)) + -a. ^^0^ + ma. ((^^)) + §• (62) 

Multiplying (62) by — li and adding the result to the right hand side of (56) yield 

E^uj- -pu" + a^payl^^ - miia. / \ \ - nVE^ + _ ug + ^(M). (63) 

The last term is considered in Appendix 2. Substituting the time derivatives J and p by 
means of (53) and (51) one obtains after cumbersome calculations: 

-E + Vfc \^pi?Uk + Ui (Ilife - piiiiik + Eodik + {{ef))5ik) 

du ■ 
dxk 

de 



iiijoi - UkJoi } + {eif) + (f^/) + W-ug (64) 



where 



A.fc = 2^a,< a,fc = -2^a,<. (65) 
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The last three terms in (64) describe the change of the field energy, field momentum and 
the effect of external forces. They depend on the concrete type of interaction. 

If there are no external fields then the last three terms in (64) should be omitted and 
the energy and momentum fiuxes are: 

(dH \ 1 
^-Q^f) - + ^ikUk (66) 

(Oe Oe \ I dH dH \ 
'^dpdp^/'^^Y'dp^dpi^/' ^^^^ 

The momentum fiux tensor consists of two parts, Ljj. and Pj^, which correspond to the 
contributions of the lattice and quasiparticles respectively: 

Lik = -{cTik + EoSik) + PoUiiik (68) 
P,HAU)+-(/f|^) (69) 

where 

is the quasiparticle momentum fiux tensor in the system of centre of mass while (Xikf) 
corresponds to the co- moving frame. It can be shown (Appendix 3) that the sume aik + EgSik 
corresponds (but coinsides in linear approximation only) to the stress tensor of the linear 
elasticity theory and turns into pressure for isotropic media. 

Finally, the equation of the elasticity theory for an elastic crystalline body with quasi- 
particle excitations in the absence of external fields takes the form: 

The last term in the right-hand side describes the force which appears when varying the 
quasiparticle mass current with respect to the lattice. 

Let us now consider the effect of the electromagnetic field. The Maxwell's equations are: 
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divD = Anq, divB = 0, (72) 

where je = jg + qu, q = qo — en is the charge density (go being the lattice charge) and 



Je = -^jo, e>0 (73) 

is the electron current density in the co-moving frame. 

We have now to take into account the field terms in (64) containing the densities of the 
field energy W, the field momentum g and the Lorentz force F. 

Before going on we would like to mention that the extraneous charges qo have to be 
considered as a separate system. They are not accounted in the Boltzmann equation as well 
as in the Lorentz force (40). Therefore, one needs some additional equations. In metals, this 
additional equation is the quasineutrality condition @] go = en. This is a good approximation 
also for semiconductors at low frequencies. The behaviour of the systems of electrons and 
other charges in the crystal depends on the problem considered. This is not the aim of our 
work. The presence of extraneous charges makes the whole system open and conservation 
laws (47)-(55) do not present a full system. One has to take into account both mechanical 
work and that of the induced forces. If the only current carriers are electrons, then one has 
to put q = 0. However, in order to keep the general form of Maxwell's equations we shall 
compensate the missing contribution Q to the time derivative of the field energy by a term 
w = -u (qE + ^je X . 

Therefore, the full field contribution is: 

Wf = ^F^/^ + W-xig + w. (74) 
Maxwell's eqs. (71) and (72) are written in L-system. They have to be supplemented 



^ Actually, if one is only interested in stress-tensor components then charges and current can be 
omitted (c.f. Ref. |o), §34). 

17 



by constitute relations. However, one has to keep in mind that these relations have their 
known simple form only in the co-moving frame. In that frame one has: 

D' = eE', B' = /xH', (75) 

where 

D' = D + -V X H, E' = E + -V X B (76) 

c c 

B' = B - -V X E, H' = H - -V X D (77) 

c c 

and V = u is introduced for convenience. 
It is easy to see that 

vD = vD' = evE' = evE, vB = vB' = /xvH' = /xvH. (78) 

Relations (75)-(78) are exact althogh (76) and (77) coincide in letter to the field trans- 
formations with an accuracy to f/cR 



^ The exact field transformations are 



E = 7E + (1 - 7)(E.n)n + 7- x B, B = 7B + (1 - 7)(B.n)n - 7- x E 

c c 

D = 7D + (1 - 7)(D.n)n + 7- x H, H = 7H + (1 - 7)(H.n)n - 7- x D 

c c 



where 



-1/2 



These fields obviously satisfy relations (78). Then from the constitute relations in the co-moving 
frame D = eE, B = ^uH one obtains immediately relations (75)-(77) in the form: 

D-F-vxH = e|E + -vxBl, B--vxE = ^|h--vxD 
c I c 1 c I c 
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When taking the time derivative of the field energy W in (74) one has to keep in mind 
that the permeabilities fj, and e in a nonstationary deformed media are functions of space 
and time. For example, the variation of the electrical part of the field energy in the lattice 
frame is: 

5W'e = ^E'5B' = |E'=^5e + eE'^E'j = —Se + S—. 
^ 47r 47r>- J Stt Stt 

Therefore, the full variation of W in time can be written in the form: 

W = — (ED + HB) = 

= fE'- ^ X bVd'- X X hU^ fH' + ^ X dUb'- X X E 

47rV c /V c J An \ c J \ c 

+ — e + — /i + v(g + G) + 



dt Stt Stt Stt 

ae(vxE')' + /x(vxH')2 (vxE')'. (v x H')^ s / 



where 



E X H ^ D X B 

Prom here on we shall restrict our consideration within an accuracy to v/c (neglecting 
terms 0{v^/c^). Then 

The quantities e and are functions of the metrical tensor g°'l^. So their time derivatives 
e and fi can be treated in the same way as the derivative of Eq (c.f. (57)-(61)). This yields: 

Ov ■ Ov ■ 

e = -^ik-K^ - vVe, /i = -i^ik-K-^ - vV/x, (82) 

OX); OX); 

where 

Substituting (82) and (83) in (81) and making use of the Poynting theoreme one obtains 



19 



/F'2 tt/2 \ Q 

W-vg= -div S' - j',E' - ^e,k + ^l^ik ^ ' 



TT Sir J oxk 
_Ve + _VMj+vG. (84) 

Neglecting terms of the order 0(f ^/c^) in (79) means that one can replace G by G'. The 
time derivative G' can be transformed using Maxwell's equations in the co-moving frame. 
This yields: 



vG' = -V (gE' + -j; X B') + ^;,V,4 + v Ve + — V/x j (85) 



where 



is the Maxwell's stress tensor in the co-moving frame. It is easy to show using eqs. (76)- (78) 
that 

V (qE' + -i'e xB')=v (qE + -je X b] = j,E - j^E'. 



J \- c 



The term which contains Lorentz force in (74) can be calculated by means of (40): 

/f— /\ - E l-e—f) - -B /— X — /\ + -li X B (-e—f\ 
\ dp / \ dp / c \dp dp / c \ dp / 

^jl^E+^xB^ ^j^E'. (87) 

The same term, but with a negative sign exists also in (84). Therefore, the total work 
related to the Lorentz force (the mechanical one and that of the electromotive forces) equals 
zero as it should be. The terms related to the extraneous charges cancel for the same reason. 
Finally, 

dv ■ 

Wf = -div S' + WkVit',, - (88) 

where 

= ^ Le'.E', + ^(e,fe - e5ik) + liH[H'^ + ^(//,fe - 1x5 A (89) 
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Hence, one has to add the term 



Q{ = Sl-Vkt[, = S,-v,W (90) 

to the energy flux density in (66), as well as the term — to the momentum flux tensor 
Uik in (67). 

The elasticity theory equation takes then the form: 

^(p^t.) = _ ^ + ^ + !!^^ _ ^ (91) 

dt * dxk dxk dxk e dt dt ' 

The current and the electromagnetic stress tensor T' are written in the lattice frame. 

Note, that the term with the electrical current represents, in fact, the electron mass- flow 

(the momentum, associated with the current). That part of the electron mass current, which 

moves together with the latice, is included in p in the left-hand side of the equation. The 

electromagnetic momentum flux tensor in L-system has the form: 

1 E"^ 
Tik = — {eEiEk + liHiHk) + (e^^ - e5ik) — + {pik - Mfc)-^ h ViGt + VkGi - ViQk - VkQi 

(92) 

Obviously, for small velocities one can replace T/^ by the corresponding tensor Tj^ in L- 
system. 

If electrons in a good conductor (metal) are considered, then the quasineutrality condition 
holds and the displacement current as well as the fleld momenta g and G have to be put 
zero. Atually, the two terms in the right-hand side of the Ampere's law can be considered 

as expansion with respect to the electric fleld frequency: j H — ~ crE + -— E. In metals 

47r ot 4:71 

a ^ c<je/(47r), and the dissplacement current can be neglected. As a result, one has 

j^rnetal ^ ^^^^^ ^ (^^^ _ ^ (93) 

47r OTT 

This tensor contains an additional term -- — /ij^ compared to our previous works 

OTT 

in which the magnetic permeability has been taken constant. As shown in Appendix 3 in 
case of a noncrystalline body (e.g. fluid) the quantities and eik have to be replaced by 



p— 6ik and p— 6ik, correspondingly (comp. Ref. ||T0|, §56). 

21 



V. CONCLUSION 



In this work we have considered dynamics and kinetics of charged quasiparticles with 
arbitrary dispersion relations in deformable crystalhne structures. We have chosen the most 
general case of time- varying deformations when the quasimomentum is not a good quantum 
number, and energy does not coinside to the Hamiltonian. We have derived a full selfcon- 
sistent set of equations which consists of the nonlinear elasticity theory equation (91) the 
kinetic equation (41) or (44) and Maxwell's equations supplemented by constitute relations. 
The kinetic equation is valid for the whole Brillouin zone and contains a new term responsi- 
ble for some inertial effects. The elasticity theory equation is derived from the conservation 
laws written in the most general form. The only approximation is that the electromagnetic 
field transformations are taken with an accuracy to 0(f^/c^). Any higher accuracy for the 
solid-state theory now would be pointless. In such a way the theory presented is exact in the 
frame of the quasiparticle approach. It can be used for any material (metal, semiconductor, 
quantum crystal, low-dimensional structures etc.) with linear relations between electromag- 
netic fields. It is easy to write the corresponding set of equations for more than one type of 
quasiparticles (e.g. electrons and holes). 

In case of electrons in metals the results obtained fit well to those of our previous works 

EE- 
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APPENDIX 1 



The evolution equation (9) for the primitive vectors can be deduced from the following 
consideration (we shall omit the subscript a for conveninence). The lattice vector a(r,t) at 
instant t is defined by the two lattice sites ri(t) and r2(t): a = r2 — ri. After a time interval 
bt the lattice vector changes to a'(r', t + = v'<^ — The new positions of the lattice sites 
are obviously r'j^ = ri + v(ri)5t, and x'^ = r2 + v(r2)5t where the velocity of the lattice site 
r in the moment t is denoted by v(r) = ii(r). It follows from the last two equations that 

a' - a = (v(r2) - v(ri))5t. 

Taking into account that v(r2) = v(ri + a) = v(ri) + (aV)v, and a' = a(r + V(5t, t + = 
a + a + (vV)a one obtains 

a+ (vV)a- (aV)v = 0. (Al.l) 

This equation coinsides with (9) in the Notation. It conserves automatically the lattice 
vectors lines. In fact, the a- vector line conservation condition consists in collinearity of a 
and the left-hand side of (Al.l), i.e. in 

[a+ (vV)a- (aV)v] X a = 0. (A1.2) 

Hence, eq. (Al.l) describes deformations which do not break or cross crystalline lines with 
equal a. This means that in a crystal lattice free of dislocations the three functions A^"(r, t) 
are single- valued, and equation (Al.l) desdribes completely the evolution of the lattice 
configuration. 

The evolution equation for the reciprocal lattice vectors a° can be obtained from (Al.l) 
and the relation 

a" = (A1.3) 

The later follows directly from the expression for the physically infinitisimal differential 
coordinates at given instant dv = aadN"' and relations (7). Multiplying (Al.l) by a^ yields 

a'^dai + a'^{uV)aai - a^a^s^ giii = (^1-4) 
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Taking into account relations (7) one has 

a,, (d^ + (uV)a^) + Vkiii = (Al.b) 

Multiplying this equation with af yields 

a" + V(a"u) = {A1.6) 

In obtaining (Al.6) we have taken into account that V^a^ — V^af = as a consequence of 
(A1.3). 

Substituting (A 1.3) into (Al.6) gives 

V (tV" + a"u) = 

and hence, 

u = -a«iV", TV" = -ua". (Al.7) 

It follows from (Al.3) and (Al.7) that 

dN" ON"' 

dN'' ^ ^—dv + ^—dt ^ a^dr - a.'^udt (A1.8) 
or at 

This expression coincides with that given in Notation. It can be written also in the form 
used in the text: 

dr = 3iadN" + udt 
APPENDIX 2 

To obtain the time derivarive of the quasiparticle energy density one can use the kinetic 
equation (44) for (f{K,r,t). Since the fluxes through the Brillouin zone boundaries in this 
notation equal zero it can be written in its integrated by parts form 

, if. 9e\ ] f_ dii de , . ^ d f , ^ 

ip+dw 1 1 u + aa^^ I (fj+aa | Ve + m— - m-^ap x curl u - F | ^-q^ = I^- (^2.1) 

We have to transform the expresssion 

^^{ef) = {{ev)) + {{e<p)). {A2.2) 
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As the derivative with respect to t is taken at constant k one can use the same procedure, 
as when obtaining (61). This yields 



+ where 



Xap = 2^ = iA2A) 



Multiplying (A2.3) by cp and (A2.1) by e one obtains after substituting into (A2.2): 
lief) ^ ((J,)) - div + a„, ((^^^))) + ^ ((f|^^)) 

= ,//)-.v(.,/).(.|/)).(F*/)-™a(|/) 

where the rules (45) are used to replace the double brakets by single ones. This expression 
is used when obtaining (64). 

We would like here to pointed out how easily this formula has been obtained. For 
comparisson, the expression which corresponds to (A2.2) in variables p, r, t has the form: 

f de\ ^ pdiii ..de ( de .. \ diii 
-Xaptti a%- mu- u(Ve)p + p,- muk 



k 



This expression is both cumbersom and nonperiodic which creates additional difficulties. 

APPENDIX 3 

The lattice contribution Lj^ to the momentum flux density tensor Ilik is given by (68). 
The term in the brakets can be written in the form: 

aik + EoSik = ((7a/3 + EoQap) a>f (A3.1) 

where Eo{g°'^) is the strain energy per unit volume. 
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In an isotropic medium the dependence of £'0 on g^"^ is reduced to a dependence on 
g = detfifa/j: 

^QEq dg 



dg dg'^f^' 
By the well known formula 

dg = -ggapdg"'^ 



one obtains easily 



and hence 



Cal3 — ■^g-Q^O'alO'Pl 



CTik = ^g-Q^aaiafiidi ctk = ^9-Q^(>ik {A3.2) 

Taking into account that g is the squared volume of a unit cell {g — V'^) and that V — 1/p 
one obtains 

CTik = ~P~Q^^^^ (A3.3) 
Hence, it is seen, that in an isotropic case 

2«>f5^ - ^^^P^-p- (^3.4) 

Let now Eq and s be the internal energy and entropy per unit mass (^0 = EqV) . 
Making use of the thermodynamic relation 

dE^ = Tds - PdV + ndN {A3.5) 

one can define the pressure P as 

It follows from (A3. 3) and (A3. 6) that in isotropic medium 

aik + EoSik = -PSik {A3.7) 

It is supposed in our consideration that the only contribution to the entropy is due 

to quasiparticles and this contribution comes from the kinetic equation. Therefore, the 
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derivatives of Eq with respect to the metrical tensor components g"'^ are assumed as taken 
at constant entropy. 

An alternative approach can be based on the free energy thermodynamic potential per 
unit volume F(T, P, N) . In that case one obtains 
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